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Abstract 

We have investigated the effects of permeable walls, modeled by linear acous¬ 
tic impedance with zero reactance, on compressible channel flow via linear stabil¬ 
ity analysis (LSA). Base flow profiles are taken from impermeable isothermal-wall 
laminar and turbulent channel flow simulations at bulk Reynolds number, Rei,= 6900 
and Mach numbers, = 0.2,0.5,0.85. For a sufficiently high value of permeability. 
Two dominant modes are made unstable: a bulk pressure mode, causing symmetric 
expulsion and suction of mass from the porous walls (Mode 0); a standing-wave¬ 
like mode, with a pressure node at the centerline (Mode I). In the case of turbulent 
mean flow profiles, both modes generate additional Reynolds shear stresses aug¬ 
menting the (base) turbulent ones, but concentrated in the viscous sublayer region; 
the trajectories of the two modes in the complex phase velocity space follow each 
other closely for values of wall permeability spanning two orders of magnitude, 
suggesting their coexistence. The transition from subcritical to supercritical perme¬ 
ability does not alter the structure of the two modes for the range of wavenumbers 
investigated, suggesting that wall permeability simply accentuates pre-existing oth¬ 
erwise stable modes. Results from the present investigation will inform the design 
of new compressible turbulent boundary layer control strategies via assigned wall- 
impedance. 


1 Introduction 


In the present work we investigate the effects of porous walls modeled with a sim¬ 
ple Darcy-like boundary condition derived as a particular case of linear acoustic 
impedance boundary condition (IBC) with zero reactance. The latter can be directly 
written in the time domain as 
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p' = Rv'n ( 1 ) 

where R is the impedance resistance, and p' and vj, are the fluctuating pressure and 
wall-normal component of the velocity at the boundary (positive if directed away 
from the fluid side), normalized with the base density and speed of sound (unitary di¬ 
mensionless base impedance). The present investigation is limited to a classic Linear 
Stability Analysis (LSA) approach and is inspired by the hydro-acoustic instability 
observed by Scalo, Bodart, and Lele ll20l . who performed numerical simulations of 
compressible channel flow with a three-parameter broadband wall-impedance (H]. 
In the latter, the effective wall permeability is a broadband function of frequency 
tuned to allow maximum transpiration only at frequency of the large-scale turbulent 
energy containing eddies. Despite the lack of the frequency-dependency in ([T]i, the 
results shown in the present investigation provide a useful theoretical framework for 
the understanding of hydro-acoustic instabilities and a guide for future investiga¬ 
tions in the manipulation of compressible boundary layer turbulence. 

The effects of porous walls on shear flows has been a topic of formidable research 
effort, especially in the low-Mach-number limit. Lekoudis m performed LSA of an 
incompressible boundary layer over two types of permeable boundaries: one mod¬ 
eling a perforated surface over a large chamber (with stabilizing effects); the second 
one, modeling pores over independent cavities (with negligible effects). Jimenez et 
ai. in] performed channel flow simulations over active and passive porous walls 
observing the creation of large spanwise-coherent Kelvin-Helmholtz in the outer 
layer responsible for frictional drag increase; the near-wall turbulence production 
cycle was not found to be significantly altered E). More recently, Tilton and 
Cortelezzi ll22l investigated channel flow coated with finite-thickness homogeneous 
porous slabs. Two new unstable modes were discovered, one symmetric and one 
anti-symmetric, originating from the left branch of eigenvalue spectrum. Similar 
results have been found in the present investigation, which also supports the find¬ 
ings of Scalo et al. ll20l who observed a complete reorganization of the near-wall 
turbulence, following the application of the tuned wall-impedance. 

One of the earliest works in linear stability analysis of compressible flow has 
been carried out by Malik na, who investigated high-speed flat plate boundary 
layers up to freestream Mach numbers of 10. Duck et al 0 performed a similar in¬ 
vestigation in the case of viscous compressible Couette flow, revealing fundamental 
differences with respect to unbounded cases. Hu and Zhong iiniini investigated 
supersonic viscous Couette flow at finite Reynolds numbers finding two unstable 
acoustic inviscid modes when a region of locally supersonic flow, relative to the 
phase speed of the instability wave, is present. Their structure was deemed con¬ 
sistent with two unstable modes shown by Mack M, which are sustained by the 
acoustic interaction between the walls and the sonic line. Malik and co-workers im 
have also used modal and non-modal stability to study the effect of viscosity stratifi¬ 
cation on the stability of compressible Couette flow. Energy transfer from the mean 
flow to perturbations occurs at two locations: near the top wall and in the bulk of 
flow domain, associated with two different unstable modes. More recent efforts are 
focused on the transitional hypersonic boundary layer stability 0123 |24] |5] 1^ . 
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where porous walls are used as means to delay transition via acoustic energy ab¬ 
sorption (allowed by a frequency-selective permeability) 

The present work aims to inform future control strategies of fully-developed 
wall-bounded compressible boundary layer turbulence over porous walls. Compar¬ 
ison against previous LSA studies on channel and Couette flow is first discussed. 
Impermeable-wall laminar and turbulent channel flow data used as a base flow in 
the present study are then presented. Finally, the structure of the unstable modes 
triggered by the porous walls is discussed along with their predicted effect on the 
structure of near-wall turbulence. 


2 Problem Formulation 

The full nonlinear governing equations are first introduced, followed by their lin¬ 
earized counterpart, used for the present LSA. All quantities reported hereafter, in¬ 
cluding results from previous works are normalized with the channel’s half-width 
(or total height in the case of Couette flow), the bulk density (constant for channel 
flow), and the speed of sound, temperature and dynamic viscosity at the wall. 


2.1 Governing Equations 

The dimensionless governing equations for conservation of mass, momentum and 
total energy are, respectively. 



( 2 ) 


(3) 


(4) 


where x\,X 2 and .*3 (alternatively x, y, and z) are, respectively, the streamwise, wall- 
normal and spanwise coordinates, «,■ the velocity components in those directions. 
Thermodynamic pressure, density and temperature are related by the equation of 
state p = j^^pT, where 7 is the ratio of specific heats. E is the total energy per 
unit mass and Rca the Reynolds number based on the reference length, (either the 
channel’s half-width or wall-to-wall distance for Couette flow), and speed of sound 
at the wall temperature, which is related to the bulk Mach and Reynolds numbers 
via Ret, = MjjRea- The viscous stress tensor and conductive heat fluxes are 
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Tij = 2ll 



1 dllfn ^ 

3 dxm 


1 II dT 
7 — 1 Pr dxj 


(5) 


where Sij is the strain-rate tensor, li the dynamic viscosity, respectively given by 
Swjj = 5 -1- , li = T’' where n = 0.75 is the viscosity exponent and Pr is 

the Prandtl number. 

The complete set of governing equations (l2]i, Q, (|4]i are only solved for com¬ 
pressible turbulent channel flow with impermeable, isothermal walls in discretized 
form on a Cartesian domain and solved with the high-order structured code Hy¬ 
brid, originally developed by Johan Larsson for numerical investigation of the fun¬ 
damental canonical shock-turbulence interaction problem OKBl. The code has 
been scaled up to 1.97 million cores in recent remarkable computational efforts 
by Bermejo-Moreno et al. El. A solution-adaptive strategy is employed to blend 
high-order central polynomial and WENO schemes in presence of strong flow gra¬ 
dient lUIIll. A forth order discretization both in space and time has been adopted 
in the present work. Results from the high-fidelity simulations are used to generate 
the base flow (see section l4Tt for the present linear stability analysis. 


2.2 Linearized Equations 

Decomposing a generic instantaneous quantity, a{x,y,z) into a base state, ,4(y), 
and a two-dimensional fluctuation, a'{x,y,t), assuming ideal gas and retaining only 
the first order fluctuations yields the following equations for conservation of mass, 
streamwise and wall-normal momentum, and energy ifTSll : 
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where Ij = j — 2/3, M = T", 1C = CpM./Pr, where M, 1C, T and U are the base 
dynamic viscosity, conductivity, temperature and streamwise velocity. In the deriva¬ 
tion of equations (| 6 ]l, (|7]i, ®, and (|9|l it was assumed that 7^ = 1 /y. The base density 
is derived from the base pressure and temperature via the equation of state with gas 
constant equal to I/ 7 , consistently with the normalization adopted. The effect of 
eddy viscosity on the fluctuations is neglected. 

Assuming harmonic two-dimensional fluctuations 




( 10 ) 


where a is the wavenumber in the streamwise direction, and (O — ac and c = Cr + i Ci 
are the complex frequency and complex wave (or phase) speed. Applying a Gauss- 
Lobatto-Chebyshev discretization along the wall-normal direction yields the gener¬ 
alized eigenvalue problem 

A‘F = ®B'F (11) 

with 'P = |u;v;p;T} where the generic column vector a is the collection of the 
discretized complex amplitudes of the generic fluctuating quantity, a'. Four condi¬ 
tions at each boundary are needed to close the system (HUi: no slip conditions for 
the streamwise velocity component, m = 0 ; isothermal conditions for temperature 
fluctuations, T = 0 ; impedance boundary conditions p = v" for y ± 1 (upper and 
lower wall, respectively); Non-homogeneous Neumann conditions for wall-normal 
derivative of pressure derived by applying ® at the boundary. 

Including the aforementioned boundary conditions on the left hand side of 
results in a singular B matrix. Therefore, eigenvalues obtained by directly solving 
the general eigenvalue problem in this form yields eigenfunctions contaminated by 
numerical noise, despite eigenvalues still being accurate. This problem can be solved 
by recasting (dUl into 

CO 'W = A ^RW (12) 

and solving the eigenvalue problem directly for the matrix C = A^^B Q. While 
this technique considerably improves the quality of the eigenfuentions, results for 
high wave-numbers (a > 2 ) still exhibit numerical issues such as top-down anti¬ 
symmetry, which has limited the current investigation to a < 2. 

The base velocity and temperature are chosen as U — u and T — T where () 
indicates Reynolds-averaged quantities, which are obtained from companion high- 
fidelity Navier-Stokes simulations (see section I tTI) . 















6 


Iman Rahbari*’* and Carlo Scalo* 


In the following, validation against LSA studies available in literature is shown 
(section|2l and new results from the present LSA analysis are discussed (sectionlDi. 


3 Comparison Against Previous LSA 

Validation is first carried out against LSA results for compressible Couette flow by 
Hu and Zhong ifTOl fTTl . where iso-thermal conditions, 7L = 1, and the tangential 
velocity, U„ are imposed at the top wall, and no-slip and adiabatic conditions at the 
bottom wall. Reynolds number and Mach number based on the topwall velocity are 
Rcoo and Moo- All quantities are normalized with speed of sound based on top wall 
temperature, the bulk density and the total wall-to-wall distance. To enable the direct 
comparison with Hu and Zhong mniiiii, the dynamic viscosity have been changed 
to: 


= C = 0.5 (13) 

and Pr = 0.72 (only for the sake of this comparison). Excellent agreement is found 
as shown by the eigenvalue spectra and eigenfunctions (figures [T] and |2]i. A grid 
convergence study of the eigenvalues (table[T]i also includes comparisons with more 
recent LSA work on the same flow by Weder ll2^ . 



Cr/Moo T, U 



Fig. 1; Comparison of complex phase velocity spectrum for compressible Couette 
flow at Mco = 2, Reoo = 2 x 10^, a = I using V=100 grid points (o) with results 

from Hu and Zhong HD (V). Real part (—) and imaginary part (-) of velocity, 

u , temperature, T , and pressure, p , eigenfunctions of the most unstable mode 
at M«, = 5, Reoo = 5 x 10^, a = 3 with (a); Re {f}, (b): Re{u\, (c): /m {f}, (d); 
Im{u}, (aa): Re {p}, (bb): Im{p}. 


Validation has also been carried out against Friedrich and Bertolotti 0 , who per¬ 
formed a stability analysis of (impermeable isothermal wall) compressible turbulent 
channel flow at Ree=4880 based on the centerline velocity, M;, = 3, a = 1 using 
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Table 1: Comparison of complex wave speed, c, of the most unstable mode for 
compressible Couette flow with previously published results for various Rec», M„ 
and a. Grid convergence study shown for A^=100, N=2QQ and N=3QQ, respectively, 
top, middle and bottom rows for each case. Deviating digits from Hu and Zhong 
data ifTOll are underlined. Values of c have been normalized with the top-wall Mach 
number. Moo, to match published values. 


Hu and Zhong do] 

~cJmZ c~JmZ 


Weder in 

Cr/M„ Ci/Moo 


Cun'ent Study 

Cr/M^ 


Re« = 5 X 10^M.„ = 5, a = 2.1, Mode I 


+0.972869314676 -0.003456356315 
+0.972869272448 -0.003456466520 
+0.972869272450 -0.003456466522 


+0.972869 178324 -0.003456358661 
+0.972869 198693 -0.003456 318849 
+0.972869201921 -0.003456324875 


+0.972869 280518 -0.0034563 23886 
+0.972869272445 -0.003456466743 
+0.972869272^ -0.0034564665^ 


Re« = 5 X 10^M.„ = 5, a = 2.1, Mode 0 


+0.040730741952 +0.000876050503 
+0.040722854287 +0.000885530891 
+0.040722853034 +0.000885531421 


+0.0407 26205831 +0.0008 84153855 
+0.04Q72285 33Q6 +Q.Q0088553 1375 
+0.040722853219 +0.000885531398 


+0.040730 596292 +0.00087 4754264 
+0.040722854373 +0.Q0Q88553Q566 
+0.040722853032 +0.000885531421 


Re« = 2x 10 -\Mo. = 2, a = 0.1, Mode I 


+ 1.213965119859 
+ 1.213965119817 
+ 1.213965119854 


-0.011585118523 

-0.011585118448 

-0.011585118558 


+ 1.213965119851 
+ 1.2139651198^ 
+ 1.213965119851 


-0.011585118549 

-0.011585118548 

-0.011585118548 


+1.213965119851 

+1.2139651198^ 

+1.213965119852 


-0.0115851185£7 
-0.011585118547 
-0.011585118547 


Re« = 2 X = 2, a = 0.1, Mode 0 


-0.291572925106 

-0.291572925140 

-0.291572925108 


-0.013821128462 

-0.013821128536 

-0.013821128457 


-0.291572925110 

-0.291572925116 

-0.291572925112 


-0.013821128465 

-0.013821128473 

-0.0138211284^ 


-0.291572925109 

-0.2915729251^ 

-0.291572925109 


-0.013821128464 

-0.013821128464 

-0.0138211284^ 




Fig. 2: Comparison of complex phase velocity spectrum for laminar compressible 
Couette flow at Moo = 5, Reoo = 1 x 10^ and a = 2.5 using A^=100 grid points (▼) 
against Hu and Zhong ifTOll . Real part (—) and imaginary part (— • —) of velocity, 
M, temperature, T, and pressure, p, eigenfunctions of the most unstable mode (indi¬ 
cated with downward arrow) with (a): Re {f}, (b); Re{ii}, (c); Im{T }, (d): Im{u}, 
(aa): Re{p}, (bb): Im{p}. 
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(bb) 




(aa) 


- 0 . 2 - 0.1 0 0.1 0.2 
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Fig. 3: Complex phase velocity spectrum for compressible channel flow at = 0.5, 
Rec = 4880 and a = 1 using N = 165 grid points (o) compared with data from 
Friedrich and Bertolotti lH (T). Real part (—) and imaginary part (-) of eigen¬ 

functions of the most unstable mode indicated with a downward arrow, (a); Re {f}, 
(b): Re{ii}, (c): Im{f }, (d): Im{u\, (aa): Re{p}, (bb): Im{p} 


Parabolized Stability Equation (PSE) with a multi-domain spectral discretization 
method. The base flow was taken from the Reynolds-averaged velocity and temper¬ 
ature profiles in Coleman et al. IH. The eigenvalue spectrum is compared in figure 
[3] where eigenfunctions of the most unstable mode (not originally plotted in their 
paper) are also included. Despite the fundamental difference between our (clas¬ 
sic) ESA approach and their PSE, the agreement is still acceptable. Eriedrich and 
Bertolotti ||8l did not explicitly mention what law for viscosity has been used, there¬ 
fore the results in figure [3] have been reproduced assuming ji = T", with n = 0.7, 
consistent with the simulations by Coleman et al. a. 


4 Compressible Channel Flow with Porous Walls 


In the following results from impermeable-wall compressible turbulent channel flow 
high-fidelity simulations are first discussed (section ED. The mean velocity and 
temperature profiles extracted from the latter are used for the present linear stability 
analysis with porous walls (section ED . 


4.1 Impermeable-Wall Turbulent Channel Flow 

The base velocity, density and temperature profiles for the present ESA investiga¬ 
tions have been taken from laminar and turbulent numerical simulations of fully 
compressible impermeable isothermal-wall channel flow. 
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The computational domain considered for the turbulent simulations in this study 
is X Ly X L, = % X 2 X A which is discretized with a number of control vol¬ 
umes Nx X Ny X = 256 x 128 x 192 resulting in a quasi-DNS resolution of 
Ax+ ~ 12.27 - 13.22, 4z+ - 8.17-8.81, =0.45-0.47, over the range 

of bulk Mach numbers investigated. The superscript + indicates classic wall-units, 
5v — u-cPwllln; where Hr = \JIwIPw is the friction velocity and p^^. and are 
the density and dynamic viscosity evaluated at the wall and the wall-shear stress is 
Tw = Pn. du/dy\„ where () indicates Reynolds averaging. 




Fig. 4: Profiles of turbulent mean velocity following Van Driest transform with ref¬ 
erence law of the wall and log-law k = 0.41 and C = 5.2 (• • •) (a) and Reynolds 
stresses normalized by the average wall shear stress, T„, (b) 


Mean velocity profiles in figure (|4]a) are plotted using Van Driest transform fol¬ 
lowing 

/■C/+ / 0 \ V2 I 

V^D= — dU+ = -\ny++C (14) 

Jo \PwJ K 

In this set of simulations, friction Reynolds number, Re^ = Urdjv, varies from 
385.6 to 417.4. Collapse of both mean velocity and Reynolds Stress profiles is ob¬ 
tained, showing that the state of near-wall turbulence has been successfully kept 
constant while varying the bulk Mach number. 

The laminar base flow has been generated with the same code with Ny = 192 
control volumes in the wall-normal direction, for values of Ret and Mi, matching 
the turbulent simulations, and are shown in flgure|2 
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Fig. 5; Profiles of laminar velocity, u{y), temperature, T{y), and density p(y) and 
Reynolds-averaged turbulent velocity, u{y), temperature, T(y), and density p(y) at 
Ret = 6900 andMt = 0.2, 0.5, and 0.85. Symbols show laminar results calculated at 
Ret = 6900 and Mt = 0.85 obtained separately by discretizing the governing equa¬ 
tions retaining only the streamwise velocity component and wall-normal gradients. 


4.2 Linear Stability Analysis With Porous Walls 

In the following we explore the effects of porous walls on the linear stability of the 
two-dimensional perturbation (fTOl i using the impermeable channel flow results gen¬ 
erated in section 14.11 as the base flow. Since turbulent eddy viscosity is neglected 
in (|2l), (O, and (|9]l, the nature of the instability between laminar and turbulent cases 
only differs due to the base flow. In particular, the profile of the wall-normal gradient 
of the mean streamwise velocity plays a fundamental role in shaping the unstable 
modes triggered by the wall permeability and determining the location of the energy 
production regions. In the following, the trajectory of the eigenvalues of the domi¬ 
nant unstable modes is analyzed, together with their induced excess Reynolds shear 
stresses. 


4.2.1 Unstable Modes Eigenvalues Trajectory 

A parametric study has been performed by varying the impedance resistance, R, 
from (zero permeability) to R = 0.01 (high permeability). Two modes. Mode 

0 and Mode I, are made unstable (c, > 0) for a sufficiently low value of resistance 
R < Ret (figure |6]l for both laminar and turbulent base flows. In all cases, the phase 
speed of the unstable modes increases monotonically with increasing the permeabil¬ 
ity, exceeding the bulk velocity. 

In the laminar case, for a = 1, the two modes originate from different regions of 
the eigenvalue spectrum. Mode I at R oo exhibits a very slow decay rate, requiring 
low values of permeability (high value of resistance Rcr = 1) to become unstable. For 
near-impermeable conditions. Mode 0 is instead initially located close to the junc¬ 
tion of Y-shape spectrum (phase velocity approximately equal to the bulk velocity) 
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Mfc = 0.2 Mb = 0.5 Mb = 0.85 






Cr — Mi, 



Cr - Mft 


Fig. 6: Trajectory of Mode 0 and Mode I for a=l in the complex phase velocity 
spectrum of compressible channel flow with impedance boundary conditions O 
using laminar (a) and turbulent (b) base flow at Ret, = 6900, M^, = 0.85 traced by 
varying the impedance resistance in the range = 10 — 0.01 for = 200 grid points. 
In increasing order of permeability, R (zero permeability): (o), R = I : D, R ~ 
0.5 :*,/? = 0.1 : V, /? = 0.05 :+,R^ 0.01 : x. 


and becomes unstable at much lower values of resistance Rcr 0.05 (higher values 
of permeability). However, the differences in their instability dynamics and spatial 
structures (discussed later) are not due to their different location in the eigenvalue 
spectrum. As a increases (not shown), the two starting points, and the corresponding 
trajectories, become closer such that for a = 2 they start from two adjacent nodes at 
the middle of the left branch of Y-spectrum. Decreasing a yields the opposite effect, 
with trajectories originating from opposite ends of the Y spectrum for a = 0.1. 

The impermeable-wall spectrum for a = 1 does not exhibit a clear Y-shape for 
the turbulent case. The trajectories of the modes follow very similar paths in the 
range of a = 0.1 —2 (not shown) tested and share very similar values of critical 
impedance resistance (or wall permeability), suggesting their coexistence. This be¬ 
havior is due to the flatness of the mean velocity profile, which is even more ac¬ 
centuated in the case of fully developed turbulence over tuned wall impedance, as 
observed by Scalo et al. 

The growth rates of the unstable modes, in case of laminar base flow, increase 
monotonically for decreasing R in the range investigated. However, when using tur- 
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bulent base flow, the growth rate saturates and then decreases for the same range 
of R. For the turbulent case only, two new unstable modes (with a similar pairing 
behavior as Mode 0 and Mode I) originate from the central part of the spectrum. 
They are not present in the laminar base flow case, travel at roughly the bulk veloc¬ 
ity, (cr Mj,) and become unstable, for very low values of resistance Rcr 0.05. A 
detailed study of them is deferred to future work. 

The transition from subcritical to supercritical permeability does not qualitatively 
alter the structure of Mode 0 and Mode I; both modes retain a very pronounced 
acoustic-like structure in the core of the channel (flgure|7]i, even in the impermeable 
wall limit. For the range of a and Mach numbers investigated. Mode 0 always man¬ 
ifests as a bulk pressure mode, causing symmetric expulsion and suction of mass 
from the porous walls, whereas Mode I resembles a standing wave, with a pres¬ 
sure node at the channel centerline. The only exception is Mode 0 in the case of 
laminar base flow profiles at subcritical permeabilities, which does not exhibit an 
appreciable structural coherence. 

Increasing the wavenumber a results in an increase of relative intensity of 
the wall-normal transpiration with respect to the centerline value of |v| for both 
modes (figure |2t) and a thinning of the boundary layer of the perturbation, ex¬ 
tending the inviscid acoustic-like core region. Future studies will verify if such 
effects are consistent with a reduction of the Stokes boundary layer thickness, 
5s = ^J2vl(0r = ^J2vl{acr). 

As permeability increases, the shape of Mode I deviates from a typical standing- 
wave-like structure, with two new pressure nodes appearing in the near wall region. 
An increase in the relative intensity of the wall-normal transpiration with respect to 
Vc in Mode 0 (figure [TJi) is also observed. Varying bulk Mach number, in the range 
considered, does not significantly alter the structure of both modes (figure]?};). 


4.2.2 Perturbation-Induced Reynolds Shear Stress 


The modifications to the Reynolds Shear Stress (RSS) distribution that would result 
from the application of porous walls can be qualitatively predicted by analyzing the 
normalized, perturbation-induced excess shear stresses 


^12 


Re{ur} 

^^^Re{uv*'\dy 


(15) 


In the case of turbulent mean flow profiles, both modes generate an anti-symmetric 
RSS distribution concentrated in base flow’s viscous sublayer (figure|4}r) with neg¬ 
ligible effects in the turbulent core (—0.9 < y < 0.9). 

When employing a laminar base flow, normalized RSS for Mode 0 using laminar 
base flow do not change considerably as permeability varies. However, the corre¬ 
sponding eigenmodes are structurally incoherent for subcritical permeabilities, as 
also pointed out in section 14.2.11 Results for Mode I show that its RSS peak ap¬ 
proaches the wall for increasing the permeability (decreasing R value). With a tur- 
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Fig. 7: Comparison between Mode 0 (left) and Mode I (right) showing pressure, 
p, (black) and wall-normal velocity, v (grey) eigenfunctions with base flow taken 
from turbulent mean profiles at Re^ = 6900. Resistance R ~ 0.01 and M* = 0.85, 

for wavenumbers a = 0.1 (—), a — 0.5 (—), and a = 1.0 (-) (a); Wavenumber 

a = 1 andMfe = 0.85, for resistances R = 1.0 (- ),R = 0.1 (—), andR = 0.01 (— 

) (b). Wavenumber a = 1 andR = 0.01, for bulk Mach numbersM^ = 0.2 (— ),Mi,= 

0.5 (-), andM^ = 0.85 (-) (c). In all Mode 0 plots, all pressure eigenfunctions 

have been normalized with their value at the centerline, pc', the resulting phase shift 
has been applied to the eigenfunctions of wall-normal velocity. In all Mode I plots, 
the vertical velocity eigenfunctions have been normalized with their value at the 
centerline, Vc, the pressure eigenfuctions’ phase has been shifted accordingly and 
then further changed by 90° to highlight the standing-wave-like structure in the 
channel core. 
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bulent base flow, the RSS shape for both modes exhibits very similar features. De¬ 
creasing the R value brings the RSS peak closer to the wall with normalized values 
evanescent in the channel core (in contrast with the laminar base flow results). It 
should be noted that increasing the permeability (lowering value of R) results in a 
higher growth rate of the disturbance, and will therefore yield a more intense aug¬ 
mentation of the total turbulent RSS (figure|8ll. 


Mode 0 



Ri2 



0 1 2 3_4 5 6 7 

Ri2 




Fig. 8: Near-wall profiles of normalized perturbation-induced RSS, Rn (ffSl l for 
Mode 0 (left) and Mode I (right) predicted by LSA for N =200 grid points, for 
different values of resistance, R, and a = 1 with base flow taken from laminar (a) 
and turbulent (b) channel flow calculations at Re* = 6900, M* = 0.85. Arrows in¬ 
dicate increasing values of permeability (decreasing value of resistance) with R = 1 

(o), R = 0.5 (••), R = 0.1 (-), R = 0.05 (- - -), R = 0.01 (—). Results for R=1 

and R=0.5 are omitted for Mode 0 with a laminar base profile, since corresponding 
eigenmodes do not exhibit a coherent structure. Turbulent RSS from the imperme¬ 
able wall calculations (figure IHi are reported here (in arbitrary units) with thick 
dashed grey lines. 


The predicted shape of the perturbation-induced RSS profiles is consistent with 
the structure of the Kelvin-Helmholtz instability observed by Scalo et al ll20l . 
sustained by the wall-normal acoustic resonance generated by the assigned wall 
impedance. 













Linear Stability Analysis of Compressible Channel Flow with Porous Walls 


15 


5 Discussion and Future Work 

In the present paper. Linear Stability Analysis (LSA) is employed to predict the ef¬ 
fects of wall-permeability on a laminar and turbulent compressible channel flow up 
to transonic bulk Mach numbers, for a fixed bulk Reynolds number of Ret = 6900. 
Two-dimensional streamwise-periodic disturbances are assumed and porous walls 
are modeled by linear acoustic impedance boundary conditions (IBC) with zero re¬ 
actance. This results in a purely real impedance, defined uniquely by its resistance, 
R, inversely proportional to the wall-permeability. The resistance has been varied in 
the range, R=\— 0.01. 

Increasing the permeability, or equivalently the inverse of the resistance R^^, 
destabilizes the flow by triggering two modes; one characterized by bulk (quasi¬ 
uniform) pressure oscillations in the core, with gradients of fluctuating quantities 
concentrated in the near wall viscous region (Mode 0); the other resembling an 
acoustic standing wave resonating in the channel’s core, with similar near-wall be¬ 
havior (Mode I). Mode 0 is characterized by symmetric u', T', and p'— eigenmodes 
and anti-symmetric v'. The latter is associated with an anti-symmetric, periodic, 
expulsion and suction of mass from the boundaries driven by the bulk pressure fluc¬ 
tuations in the core. On the other hand. Mode I exhibits anti-symmetric u', T', and 
p' profiles with symmetric v' distribution. For both modes, increasing the wall per¬ 
meability causes the perturbation-induced Reynolds shear-stress peak to shift closer 
to the porous walls, for both laminar and turbulent base flows. 

For very high values of permeability (low values of resistance R < 0.01) and 
only in the case of a turbulent base flow, a saturation of the growth rate is observed. 
Values of R < 0.1 may, however, be unrealistic (see measurements of the acoustic 
response of perforated panels ll2Tll ') and are also expected to yield nonlinear effects 
that are not captured by a linear IBC formulation. 

Although the current analysis is limited to purely real acoustic impedances (i.e. 
classic Darcy-like formulation for porous walls), the results obtained qualitatively 
confirm the structure of the hydro-acoustic instability observed by Scalo et al. IJOl. 
In their case, large spanwise-coherent Kelvin-Helmholtz rollers are found to be con¬ 
fined near the wall where Reynolds shear stresses are significantly augmented by 
the instability, while the outer layer remained unaffected. This result is in qualita¬ 
tive agreement with the perturbation-induced Reynolds shear stress distribution pre¬ 
dicted by the current analysis and in contrast with the findings by Jimenez et al. na, 
valid for the incompressible limit, where Kelvin-Helmholtz rollers were observed in 
the outer layer and the near wall turbulence was unaltered. A more accurate compar¬ 
ison between our LSA and the results by Scalo et al. ll20ll will be obtained in future 
studies, which will incorporate a broadband reactance (the impedance’s imaginary 
part) in the eigenvalue problem. 

Eigenvalue tracking shows that the modes made unstable by permeability pre¬ 
serve the characteristic structure they had when stable. Only in the laminar case. 
Mode I is at the verge of instability for zero permeability. As the wavenumber a is 
increased, trajectory of the unstable modes become close to each other, and near¬ 
wall gradients of fluctuating quantities intensify. 
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Future work will focus on Direct Numerical Simulation of channel flow with 
general impedance boundary conditions and new companion LSA efforts aimed at 
finding the optimal impedance distribution for effecting flow control, with improved 
numerics and physical modeling of the boundary conditions. Investigations will be 
also carried out considering boundary layers over flat plates, to rule out confinement 
effects. 


6 From Small Data to Big Data (and Vice Versa) 

This work has been centered around the use and generation of small data, in the form 
of a classic Linear Stability Analysis, in order to acquire a general understanding 
of the response that it is expected from compressible channel flow turbulence under 
the effects of porous walls. This problem can only be truly explored in detail with 
large-scale, high-fidelity simulations, that is, big data. Simplified low-order models, 
however, prove to be essential in guiding the design and interpretation of companion 
large-scale numerical simulations, since they provide us with a rough but robust 
understanding of the physical problem under investigation at a considerably reduced 
cost. Moreover, low-order models allow for a true parametric study, spanning, for 
example, several orders of magnitude of variation of a key parameter in the problem, 
such as, in this case, the wall permeability. High-fidelity simulations, on the other 
hand, can be used to populate a few key regions of the parameter space with limited 
modeling assumptions, and used to calibrate and/or verify the results from the low- 
order modeling efforts. 

When possible, a synergistic interaction between small data and big data should 
therefore be sought. A “blind big data”- or a “blind small data”- only approach may 
suffer from severe limitations. The development, analysis and adoption of the two 
approaches should go hand in hand in order to cross-verify results, especially when 
experimental data is not available. However, with the advances in computational 
power, the definition of small and big data is changing rapidly. For example, tur¬ 
bulence simulations performed on a computational grid of = 1024^ points, can 
nowadays (2015) be considered average-size data. In spite of the gradual shift in 
the demarcation line between small and big data, we believe that the importance of 
low-order models (used here as an example of small data) will never (and should 
never) cease to retain a prominent role in the design, analysis and support of results 
from large-scale numerical investigations. 
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